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We study the effect of uncertainties in the proton-capture reaction rates of the NeNa and MgAI 
chains on nucleosynthesis due to the operation of hot bottom burning (HBB) in intermediate-mass 
asymptotic giant branch (AGB) stars. HBB nucleosynthesis is associated with the production of 
sodium, radioactive 26 Al and the heavy magnesium isotopes, and it is possibly responsible for the 
0, Na, Mg and Al abundance anomalies observed in globular cluster stars. 
We model HBB with an analytic code based on full stellar evolution models so we can quickly 
cover a large parameter space. The reaction rates are varied first individually, then all together. 
This creates a knock-on effect, where an increase of one reaction rate affects production of an 
isotope further down the reaction chain. We find the yields of 22 Ne, 23 Na and 26 Al to be the most 
susceptible to current nuclear reaction rate uncertainties. 
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1. Introduction 

The Asymptotic Giant Branch (AGB) is the final evolutionary phase of stars with a mass less 
than about 8 M Q . The core of an AGB star has exhausted its supply of hydrogen and helium by 
nuclear burning and is made of carbon, oxygen, neon and perhaps magnesium, the main helium- 
burning products. A thin helium layer separates it from the convective stellar envelope which is 
mostly hydrogen. Hydrogen burning in a shell at the base of the envelope dominates the stellar 
luminosity, which is many thousands of times that of the Sun, but is punctuated by thermal pulses 
due to ignition of the helium layer. At each pulse the third dredge-up may occur, in which helium- 
burnt material is mixed into the stellar envelope, polluting it with helium, carbon, 22 Ne and heavy 
s-process elements. The lifetime of an AGB star is dominated by rapid mass-loss which limits it to 
about a million years, after which the stellar envelope is ejected into the interstellar medium and a 
CO white dwarf remains. 

In AGB stars more massive than about 4M Q the hydrogen burning shell, at a temperature of 60 
to 100 million K, extends into the convective envelope, a process called hot-bottom burning (HBB). 
The envelope composition is directly affected by the various hydrogen-burning cycles: CNO, NeNa 
and MgAl. The effect is to convert carbon and oxygen to nitrogen, neon to sodium and magnesium 
to aluminium. The abundance anomalies seen in globular cluster stars, for example the oxygen- 
sodium anticorrelation, may be due to HBB in AGB stars []]J]. The combination of third dredge-up 
of carbon with HBB is an important source of primary nitrogen, especially at low-metallicity. 

Stellar models of hot-bottom burning rely on nuclear reaction rates that are inherently uncer- 
tain due to both experimental limits and extrapolation into the low-energy stellar regime. Our aim 
is to determine the effect of considering the rate uncertainties on the stellar abundances, and hence 
the chemical yields from AGB stars undergoing HBB. We have selected the NeNa and MgAl pro- 
ton capture reactions primarily because these are modelled directly by our synthetic AGB model. 
It is not possible to vary the CNO reactions because they contribute significantly to the stellar lu- 
minosity, so affect the stellar structure, while we are assuming that changes in the NeNa/MgAl 
reaction rates do not affect the stellar structure. 

2. Model 

Our synthetic AGB model is described in detail in [Q] but has been extended to include the 
NeNa and MgAl cycles (Izzard et al 2006, in preparation) as well as the CNO cycle. We model 
burning in the convective envelope by replacing the burn/mix/burn/mix. . . cycle with a single 
burn/mix event (see Figure [[]) where the amount of material burnt, and the time it is burnt for, is 
calibrated to detailed AGB models [0]. 

We analytically solve the differential equations for the abundances of 20,21 ,22 Ne, 23 Na, 24 ' 25 26 Mg 
and 26g,26m,27 A1 in the burning region so we can follow the surface abundances as a function of time. 
The NeNa cycle is solved by an eigenvalue method (assuming no input to or output from the cycle, 
which is justified as 23 Na(/?, y) 24 Mg is usually slow and input from 19 F is negligible compared 
to the amount of 20 Ne present), then 23 Na is decayed explicitly into 24 Mg. The MgAl cycle then 
consists only of forward reactions if we ignore 27 Al(p, a) 24 Mg, which is always justified at the 
temperature of HBB (log 10 r/K ~ 7.7 — 8.0), so each rate equation is easily solved analytically. 
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Figure 1: Schematic view of our hot-bottom burning model. The detailed stellar models show that only a 
thin layer at the base of the convective envelope is hot enough to burn. In reality, and in the detailed models, 
the thin layer is continously mixed over hundreds of convective turnover timescales, but our synthetic model 
can reproduce the effect by burning a large fraction of the stellar envelope for a short time. 



Reaction 


Lower limit 


Upper limit 


Reference 


20 Ne(/j,7) 21 Na(j8 + ) 21 Ne 


-50% 


+50% 


Iliadis et al. 2001 ^ ]=NACRE | J] 


21 Ne( j p,7) 22 Na( j 8+) 22 Ne 


-20% 


+20% 


madisetal. 2001 [ 




22 Ne(p,y) 23 Na 


-50% 


x2000 


Hale el al. 2001 [[J 


23 Na(/j, a) 20 Ne 


-30% 


+30% 


Rowland et al. 2004 Q 


23 Na(/?,y) 24 Mg 


/AO 


xlO 


Rowland et al. 2004 




24 Mg( j p,7) 25 Al( j 8+) 25 Mg 


-17% 


+20% 


Iliadis et al. 2001 [|] = Powell et al. 1999 | J 


25 Mg(/j,7) 26 AF( j 8+) 26 Mg 


-50% 


xl.5 


Iliadis etal. 2001 [|| 


26 Mg(/?,7) 27 Al 


/4 


xlO 


niadisetal. 2001|5 




26 Mg(/?,7) 27 Al 


-25% 


x3 


Iliadis etal. 200l{]] 


26 AF(/7,7) 27 Si(j8+) 27 Al 


-50% 


x600 


Iliadis etal. 200l(]] 



Table 1: Nuclear reaction rates in our network and their uncertainties. 



Beta-decays except that of 26 Al g are treated as instantaneous, we do allow for 21 Al(p, a) 24 Mg by 
explicit decay of 27 Al (the rate of this reaction is always negligible) and the chain is terminated 
at 27 Al because 27 Al(p, 7) 28 Si is very slow. The nuclear reactions, rate uncertainties in the HBB 
temperature range and references are found in Table [T[ 
We have constructed two sets of models in which: 



1 . Each reaction is varied individually, 

2. All reactions are varied at the same time. 



When comparing the yield of an isotope, which we define simply as the mass ejected as that isotope, 
we use a control model where all the reaction rates are set to their recommended values (i.e. their 
rate multipliers are all 1.0). We useM = 5,6M and Z = 0.02,0.004 to test the effect of increasing 
mass and decreasing metallicity, both of which increase the amount of hot-bottom burning. 



3. Results 

Figure |2| shows the effect on chemical yields of the NeNa and MgAl isotopes of varying the 
proton-capture nuclear reaction rates. The values shown are percentage differences from the control 
values, where positive values indicate extra production, negative values extra destruction, and zero 
means no change. 

The main effects are: 
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Figure 2: Uncertainty ranges for chemical yields as the difference from the control value as a function 
of isotope for our four sets of model. Dark colours are the maximum uncertainty ranges for varying each 
reaction individually, the light colours give the uncertainty when all the reactions are varied together. 

• Complete conversion of 22 Ne to 23 Na, or not, due to the large uncertainty in the 22 Ne(p, y) Na 
rate. The yield of 23 Na can be boosted by a factor of about 40 at Z = 0.004, in both the 5 and 
6M models. This sodium is primary; 

• A very large uncertainty in 26 Al production, ranging from x to x 2. This stems from the 
x600 range of the 26 Al(p, y) 27 Si rate. There is a corresponding knock-on effect on 27 Al; 

• Most other isotopes vary by at most about 10 — 20%. 

It is sometimes necessary to consider changing all the reaction rates together. The yields of the 
magnesium isotopes at Z = 0.004 show this clearly, as they are affected by the preceding NeNa 
cycle rate changes followed by 23 Na(>, y) 24 Mg ...{p, y) 25 Al(/3+) 25 Mg ...(/?, y) 26 Al(/3+) 26 Mg. 

4. Conclusions 

We have showed that changes in yields from hot43ottom burning stars due to nuclear reac- 
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tion rate uncertainties can be significant, particularly for Ne, Na, the magnesium isotopes and 
aluminium. The effect of reaction rate uncertainty becomes stronger as the hydrogen burning tem- 
perature increases i.e. at higher mass and/or lower metallicity. For some isotopes, especially those 
of magnesium, it is important to consider changing all the uncertain reaction rates at the same time, 
as simply changing the proton capture rate on each isotope does not give an accurate error range 
for the yield. 

It remains to be seen if, say, a higher Ne(p, y) 23 Na rate can help to solve the problem of the 
globular cluster abundance anomalies, which show excessive sodium and lower oxygen, particu- 
larly as we have not tried to change the CNO cycle rates here. It is also possible that the yields of 
primary sodium could be constrained by galactic chemical evolution models, which in turn would 
give us an upper limit for the 22 Ne(/?, y) 23 Na rate, although we will also have to examine carefully 
the role of massive stars in 23 Na production. 

Factor of two changes in the magnesium yields in our low-metallicity models may be important 
for fine-structure constant (a) variation studies because Aa/a is a steep function of ( 25 Mg+ 26 
Mg)/ 24 Mg (see Fig. 1 of [§]). 
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